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ABSTRACT 



Analyses were performed for static and dynamic buckling of 
a continuous fiber embedded in a matrix in order to determine 
the effects of interfacial debonding and fiber breakage on the 
critical buckling load and the domain of instability. A beam 
on elastic foundation model was used for the study. The study 
showed that a local interfacial debonding between a fiber and 
a surrounding matrix resulted in an increase of the wavelength 
of the buckling mode. An increase of the wave length yielded 
a decrease of the static buckling load and lowered the dynamic 
instability domain. 

In general, the effect of a partial or complete 
interfacial debonding was more significiant on the domain of 
dynamic instability than on the effects of static buckling 
load. For dynamic buckling of a fiber, a local debonding of 
size 10 to 20 percent of the fiber length had the most 
important influence on the domains of dynamic instability 
regardless of the location of debonding and the boundary 
conditions of the fiber. For static buckling, the location of 
a local debonding was critical to a free -simply supported 
fiber but not to a fiber with both ends simply supported. 
Fiber breakage also lowered the critical buckling load 
significantly. 
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I. 



INTRODUCTION 



Buckling is one of the important failure modes of 
structures subjected to compressive inplane loads. Buckling in 
a continuous fibrous composite panel may include global 
buckling of the panel as well as buckling of fibers embedded 
in a surrounding matrix. Fiber buckling is one of the major 
causes of reduced compressive strength of a composite panel. 
As a result, the compressive strength of a composite is 
generally lower than its tensile strength. 

Some of the initial studies on fiber buckling were given 
in references (Dow and Grunt fest 1960, Fried and Kaminetsky 
1964, Leventz 1964, and Rosen 1965) . These studies found that 
fiber buckling was an important cause in reducing the 
compressive strength of composites. Rosen (1965) investigated 
the fiber buckling using an energy method. He considered 
two-dimensional arrays of fibers embedded in a matrix 
material, and computed critical fiber buckling loads for two 
different modes of buckling: the shear mode and the extension 
mode. In the shear mode, two neighboring fibers buckle in 
phase. On the other hand, two neighboring fibers buckle out of 
phase in the extension mode. Of the two fiber buckling modes, 
the one which results in a lower critical buckling load 
governs the buckling mode. 
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Experimental studies were performed to investigate fiber 
buckling within a supporting matrix material (Lager and June 
1969, and Dale and Baer 1974) . The experimental results agreed 
with Rosen's predictions. For a very low fiber volume 
fraction, the extension mode predominated, while the shear 
mode predominates for composites with fiber volume fractions 
of interest in practical application. A postbuckling study of 
fibers was performed by Maewal (1981) . He found that the 
effect of initial geometric imperfections like initial 
waviness of fibers on the microbuckling stress was not 
significant . 

Experimental studies showed that compression failure in 
straight fiber laminated test specimens initiated at a free 
edge (Hahn and Williams 1984) . In order to enhance the 
understanding of the effect of the free edge on the initiation 
of the buckling process, a study of fiber buckling was 
undertaken by Wass, Babcock, and Knauss (1989) . They used a 
model of a beam supported by an elastic matrix. They concluded 
that, for low fiber volume fractions, the critical strain 
values at a free edge were lower than the predictions using 
Rosen's model. For high volume fractions, however, they found 
the beam model was not valid for the purpose of their study. 

All the studies mentioned above assumed perfect bonding 
between a fiber and its surrounding matrix. Lanir and Fung 
(1972) studied the effect of interfacial separation on fiber 
buckling. They calculated deflections of the postbuckled 
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wave-like shape of fibers. They also stated that the critical 
buckling load of a fiber was the same as that of a free column 
in case of interfacial separation between the fiber and 
matrix. This was the situation where a surrounding matrix did 
not affect the fiber buckling due to separation. However, if 
there is an interfacial separation between a fiber and a 
surrounding matrix but they are in interfacial contact without 
pressing each other before buckling, the matrix material still 
supports the fiber after the buckling occurs as shown in 
Figure 1. In a buckled position, the compressed matrix side 
(the bottom side in Figure 1) provides reactions to a fiber. 
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The present thesis computes critical buckling loads of 
fibers with partial interfacial debonding from surrounding 
matrices using the finite element method. The buckling loads 
with debonding are compared to the buckling loads without 
debonding for various fiberous materials. Different locations 
and magnitudes of partial debonding are considered to 
investigate their effects on the critical buckling load. The 
study includes not only static buckling but also dynamic 
buckling of a fiber embedded in a matrix material. 
Furthermore, the effect of fiber breakage on the fiber 
buckling is also investigated. 
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II. ANALYSIS AND MODELING 



A. NUMERICAL ANALYSIS 

A beam on elastic foundation model was used for this 
study. The Bernoulli -Navier beam theory was used to model a 
fiber, and a surrounding matrix was modeled as an elastic 
foundation. The spring constant of an elastic foundation was 
computed from the expression developed by Herrmann, Mason and 
Chan (1965) . The spring constant is given below: 



k= ( 1 ) 

(3-4u m ) K 0 (2nx/S) + (nx/b) K x (2nx/b) 

in which 

k = spring constant of the elastic foundation, 

G = shear modulus , 
v = Poisson's ratio, 
r = radius of the fiber, 

5 = wavelength of the deformed fiber (see Figure 2) , 

K 0 = zeroth-order modified Bessel function of the second kind, 
K x = first-order modified Bessel function of the second kind, 
and the subscript m denotes the matrix. The spring constant in 
Equation (1) was obtained assuming that the moment foundation 
modulus associated with the beam rotation was negligible. This 
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assumption is valid when the shear modulus of a fiber is much 
larger than that of a matrix. 



The equation of beam bending is : 




( 2 ) 



where 

m = mass of the fiber per unit length, 
w = transverse deflection of the fiber, 

E = elastic modulus, 

I = moment inertia of cross-section, 

P = axial load applied to the fiber (a positive value denotes 
a compressive force) , and 

x direction is along the fiber axis while t and f denote time 
and fiber respectively. Applying the method of weighted 
residual along with the weak formulation to equation (2) 
gives ; 
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Figure 2. Wavelength of a buckled fiber 
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In this equation, w, x indicates a partial derivative of 
w with respect to x, w is a test function. This equation 
contains second derivatives so that it is needed to assume an 
interpolation function for w that has both continuous 
deflection, w, and continuous slope, w' , between elements. One 
common choice is the cubic Hermitian polynomial presented in 
below. 



H x (r ) =2r 3 -3r 2 +l 

H 2 (r ) = (r 3 -2r 2 +r ) h (4) 

H 3 (r) =3r 2 -2r 3 
H 4 (r ) = (r 3 -r 2 ) h 

where r=x/h and h represents the length of element. Their 
derivatives with respect to x are 



H 1 / =(6r 2 -6r) / h 
H 2 =3r 2 -4r +1 
H 3 /= (6r -6r 2 ) / h 
H 4 /= 3r 2 -2r 



H 1 // =(12r-6) / h 2 
H 2 //= (6r-4) / h 
H 3 //= (6 -12r ) / h 2 
H 4 // =(6r-2) / h 



(5) 



where prime, ' , denotes a derivative 

Substituting the Hermitian shape function to Equation (3) 
after discritization of the domain into a number of finite 
elements results in: 



Jm[H] T [H]dx +j*E f I f [H // ] T [H // ]dx 

h h 





= 0 



( 6 ) 



or 
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[M] (w) + [K ei ] (w) - [K p ] (w) + [K k ] (w) = 0 . 



(7) 



The superimposed dot denotes temporal derivatives. The details 
of mart ices in Equation (7) are given below: 



w -p / 



H 1 // H 1 // h/'h" hI'h" h/'h" 

U H IJ H TJ //tT H TJ //tt // 

n 2 xij_ n 2 ii 2 ii 2 3 n 2 n. 4 

TJ ^TT H IT H-U N TT H-U H TJ ^TJ H 

ri 3 ri 2 n 3 n 2 ii 3 n 3 rl 3 rl 4 

TT //tj // TJ Wit H TJ Wtj // TJ //tj // 

n 4 n x n 4 112 ri 4 xl 3 ri 4 ri 4 



dx 
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h 3 
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-12 -6h 12 -6h 

6h 2h 2 -6h 4h 2 
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dx 
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-13h 
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22h 
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13h 
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420 
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[Kp] = 





TT ! IT ! TT ^ TJ ^TJ ^ 

nin 2 ri^ri3 ri^ri^ 


h 

P f 


Tj ^TJ ^ TJ ^TJ ^ TJ ^TJ ^ TJ ^ 

r± 2 ^i -n- 2 ri 2 rl 2 rl3 rt 2 n 4 


h 2 J 

11 0 


TJ ^ TJ ^TJ ^ TJ ^TJ ^ TJ ^TJ ^ 

113x11 11311 2 113113 H3II4 




TJ ^TJ ^ TJ ^TJ ^ TJ ^TJ ^ TJ ^TJ ^ 
n.4 ill ri4rl 2 114113 H4H4 



dx 



( 8 ) 



(9) 
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p 

30h 



36 3h -36 3h 



( 10 ) 



[M] = m 



/ 
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In Che present study, the diagonal mass matrix was also used 
in addition to the consistent mass matrix given equation (11) . 
The diagonal mass matrix is 



[M] = — 
78 



39 0 0 0 

0 h 2 0 0 

0 0 39 0 

0 0 0 h 2 



( 12 ) 



For the static buckling analysis, the mass term in 
Equation (3) drops out, and the resultant equation becomes an 
eigenvalue problem as given below: 

det [ [K ei ] + [K k ] - P CI [K p ] ] = 0 (13) 

where P cr is the critical buckling load. A compressive load 
was assumed to be positive. The elastic foundation increases 
the critical buckling load significantly compared to the case 
of a beam without any foundation support. Furthermore, mode 
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shapes of buckling become complicated due to the elastic 
foundation (see Timoshenko and Gere 1961, and Hetenyi 1974) . 
For example, the mode shape for the smallest critical buckling 
load is dependent on the geometric and material properties of 
the beam and foundation. The mode shape is unknown a priori. 

In order to determine the spring constant from Equation 
(1) , the wavelength must be known. However, the wavelength is 
unknown because the mode shape is unknown a priori. 
Therefore, an iteration process is necessary to find the 
correct wavelength of the buckled beam. The flowchart for the 
iteration procedure is given in Figure 3. 

B. NUMERICAL MODELLING 

Interfacial debonding between a fiber and the supporting 
matrix is modeled as explained below. It is assumed that there 
is no gap at the interface between the fiber and the matrix 
even if there is debonding at the interface. As a result, if 
a fiber buckles at the debonded interface, the compressed side 
of a matrix material still supports the fiber while the other 
side of the matrix does not provide any support as shown in 
Figure 1. If there is an enough friction to prevent sliding at 
a debonded interface which is located at the compressed side 
of a matrix, the interface under compression behaves just like 
a bonded interface. A fully bonded matrix along the 
circumference of a buckled fiber supports the fiber at both 
the compressed and elongated sides of the matrix, while a 
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matrix debonded along the circumference of a buckled fiber 
supports the fiber only at the compressed side of the matrix. 
Therefore, it is reasonable to assume that the matrix at a 
debonded interface provides half of the support of a perfectly 
bonded matrix if a correct wavelength is obtained at the 
debonded interface. However, the wavelength changes with the 
debonding. As a result, an iteration is required to determine 
the correct wavelength and the spring constant at a debonded 
interface . 

The fiber breakage is modeled in the following way. Once 
there is a fiber breakage, there is no continuity of the slope 
of the deflection at the broken position. However, the lateral 
deflection is assumed to be continuous. Otherwise, there is no 
axial force to balance the two broken fibers in the axial 
direction since the spring does not constrain the axial 
motion. Debonding is also assumed to occur from the broken 
position. 

For the study of dynamic instability, the axial force in 
Equation (3) is assumed to be a harmonic function. That is, 

P = p cos (0 t) (14) 

where p is the magnitude of the pulsating inplane force and 6 
is the circular frequency of the pulsation. Let the transverse 
displacement vector be expressed in terms of trigonometric 
forms; 
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( BEGIN ) 




Figure 3. An iteration procedure to find the wavelength of the 
buckled shape. 
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( 15 ) 



(w) = £ { Uj sin + y { a } cos i0t 

i*i, 3, . . 2 1 2 

where u and u are time independent vectors of the nodal 
transverse displacement. Substituting Equations (14) and (15) 
into Equation (7), and considering the first term (i.e. i=l) 
of Equation (15) as the first-order approximation yields the 
following eigenvalue matrix equation: 

det [ [K EI ] + [K K ] ± £ [K p ] - 1 = 0 (16) 

The eigenvalues of 6 obtained from Equation (16) determine the 
boundary of domain of principal instability. For the mass 
matrix in Equation (16) , both diagonal and consistent mass 
matrices are used for comparison. The transverse shear effect 
and rotatory inertia effect are neglected in this study 
because the length of the fiber under consideration is much 
larger than the cross-sectional dimension. A previous study of 
dynamic instability of layered composite plates (Kwon 1991) 
showed that the transverse shear effect was found to be 
negligible if the length of a plate was much larger than a 
plate thickness. 
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III. RESULTS AND DISCUSSIONS 



A. STATIC BUCKLING 

Under consideration in this study were fibers made of 
graphite, Kevlar, and glass embedded in epoxy, and boron 
embedded in an aluminum. Their material properties were 
different so that the effect of different materials on 
buckling was studied. The material properties are tabulated in 
Table 1, which were obtained from a textbook (Tsai and Hahn 
1980) . 

TABLE I. MATERIAL PROPERTIES OF FIBERS AND MATRICES 



MATERIAL 


ELASTIC 
MODULUS (GPa) 


DENSITY (g/cm 3) 


GRAPHITE 


230 


1.75 


KEVLAR 


120 


1.44 


GLASS 


72 


2.6 


BORON 


410 


2.6 


EPOXY 


3.9 




ALUMINUM 


70 





The first study considered a fiber of dimension L/r = 100 
embedded in an epoxy. Here L/r is the ratio of the beam length 
to radius. Both ends were assumed to be simply supported. 
Interfacial debonding was assumed to occur at the center of 



14 
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K L )| 



Figure 4. A simply supported fiber with interfacial 
debonding at the center. 

the fiber as shown in Figure 4. The effect of the interfacial 
debonding on the static buckling load is shown in Figure 5 for 
four different cases. The results revealed that the percentage 
reduction in the critical buckling load was almost independent 
of the elastic moduli of the fibers even if the absolute 
magnitude of the buckling loads varied for the three different 
fibers. The graphite fiber had the largest magnitude of 
buckling load and the glass fiber had the smallest buckling 
load because of their elastic moduli. There was approximately 
a 35 percent decrease of the static buckling load due to 
complete debonding. In the case of the boron fiber embedded in 
an aluminum matrix, there was a 42 percent decrease in the 
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Figure 5. Effect of interfacial debonding on the buckling 
load for simply supported fibers. 



static buckling load due to complete debonding. For all cases, 
it can be seen that, the reduction of the buckling load was 
significant for initial 10 to 20 percent of debonding area and 
changed very gradually after that. 

The static buckling mode shapes varied as the interfacial 
debonding increased. The shape changed from a symmetrical mode 
to an antisymmetrical mode or vice versa as the debonding 
progressed. In addition, the wavelength changed with the 
debonding. Figure 6 shows buckled shapes of the Kelvar fiber 
for different debonding sizes. The wavelength at the debonded 
area was larger than that at the bonded area. The relative 
magnitude of the buckled shape was much larger at the debonded 
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area than at the bonded area if there was a partial debonding 
through the length of the fiber. Therefore, even if there is 
waveness at the bonded area in Figure 6 (b) , it is not 
reflected in the buckled shape. The graphite fiber had the 
largest wavelength while the glass fiber had the smallest 
wavelength in an epoxy matrix. 

In order to examine the effect of fiber length on 
buckling, a fiber with ratio L/r = 200 was studied. The 
absolute magnitude of critical load decreased for a fiber of 
L/r = 200 when compared to that of a fiber of L/r = 100, but 
the percentage reduction of the critical loads caused by the 
interfacial debonding remained almost the same for the two 
fibers of different lengthes. The wavelength was larger for 
the longer fiber. The location of partial debonding was varied 
to investigate its effect on the buckling load. For a fiber 
simply supported at both ends, the location of debonding did 
not make much difference on the decrease of buckling load. An 
initial debonding at any location was the most critical factor 
in reducing the critical buckling load. 

The next study was to determine the effect of different 
boundary conditions. A fiber with one end simply supported and 
the other end free was considered. As before three different 
fibrous materials were studied. First, a partial debonding was 
assumed to start from the free end as shown in Figure 7. 
Figure 8 shows the reduction of the critical buckling load as 
a result of increased partial debonding. An initial debonding 
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(a) Perfect interfacial bonding 



(b) 50 % interfacial debonding at the cente 




(c) complete interfacial debonding 



Figure 6. Buckling modes of Kevlar fiber. 
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(a) Debonding starting from free end 



DEBONDING 

AREA 




(b) Debonding starting from simply supported end 



Figure 7. A free simply supported fiber with different 
locations of interfacial debonding. 
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Figure 8. Effect of interfacial debonding on the buckling load 
for free simply supported fibers. (Debonding starts from free 
end) 



of 10 percent decreased the buckling load significantly while 
no further debonding affected the buckling load. When a 
partial debonding started at the simply supported end as shown 
in Figure 7, there was no change in the buckling load 
untilthere was at least 10 percent debonding as shown in 
Figure 9. That is, the location was very critical in 
decreasing the buckling load. This result was contrary to that 
for a fiber with both ends simply supported. 

The last study for static buckling was to examine the 
effect of fiber breakage along with interfacial debonding. A 
fiber with simply supported ends was considered. As before 
four different fiber materials were used. After a fiber 
breakage, partial debonding was assumed to start from that 
point. (See Figure 10) Figure 11 shows a comparison of the 
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Figure 10. Fiber breakage along with interfacial debonding. 

fibers with a breakage and without it. A fiber breakage showed 
a 50 percent decrese of the critical load, independent of 
fiber breakage locations except for very close distances from 
the supports. Fiber breakage near the supports is more 
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critical. Debonding effects after a fiber breakage were very 
similar to the case without fiber breakage. An initial 10 
percent debonding decreased the buckling load significantly 
and after 30 percent no further debonding changed the buckling 
load. 




Figure 11. Comparison of critical buckling loads for a broken 
fiber and a fiber without breakage. 



B. DYNAMIC BUCKLING 

The domains of dynamic instability were plotted in Figures 
12-15 for four different cases (Graphite -epoxy, Kevlar-epoxy, 
Glass-epoxy and Boron -aluminum) . Each differnt cases has 
L/r=100 ratio with simply supported ends. The axial force and 
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circular frequency in all figures were normalized as shown 
below : 



• NORM 



P 



CR 



(17) 



0 



NORM 



6 

2 Q 



(18) 



where P CR is a static buckling load and Q is the natural 
frequency. A local debonding was assumed to begin from the 
center. An interfacial debonding lowered the instability 
domains. The results showed that the instability domains of 
fibers embedded in an epoxy matrix were lowered 35 percent at 
full debonding compared with no debonding. The instability 
domain of boron -aliminum was lowered by 42 percent for the 
same comparison. An initial 10 to 20 percent local debonding 
influenced the instability domains significantly but a further 
increase of debonding size had very little affect. Local 
debonding also widened the instability domains. Although the 
instability domains between 10 percent debonding and 100 
percent debonding were not shown in Figure 12, these 
instability domains were calculated and were located between 
the outer borders of the mentioned domains. 

After calculating the elastic moduli ratios, which was 
shown in Table 2, it was seen that, when the elastic moduli 
ratio decreased, the instability domains between 10 percent 
debonding and full debonding scattered over a smaller area. 
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TABLE II. ELASTIC MODULI RATIOS BETWEEN FIBERS AND MATRICES 



FIBER -MATRIX 


ELASTIC MODULI 
RATIO 


GRAPHITE -EPOXY 


58.97 


KEVLAR -EPOY 


30.77 


GLASS -EPOXY 


18.46 


BORON- ALUMINUM 


5.86 



When a local debonding started at both simply supported 
ends, the debonding effect was less critical than a local 
debonding starting at the center until it progressed 20 
percent debonding. For a debonding size larger than 20 percent 
of the beam, the debonding effect was not sensitive to the 
location (see Figures 16-19) . 

Fibers with free- simply supported ends were also studied 
for dynamic buckling. Figures 20 through 27 show the change of 
domains of instability caused by an interfacial debonding. 
Figures 20 through 23 represent the case of a debonding 
starting from the free end while Figures 24 through 27 
represent the case a debonding starting from the simply 
supported end. The debonding starting from the free end was 
more critical to the instability domains than that starting 
from the supported end. However, the dynamic buckling was less 
sensitive to the location of a debonding than the static 
buckling for free- simply supported fibers. The size of a 
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debonding was more important than location in influencing the 
instability domain. 
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Figure 12. Instability domains Kevlar-epoxy with free and 
simply supported ends (Dedonding started from the simply 
supported end) . 
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Figure 13. Instability domains glass-epoxy with free and 
simply supported ends (Dedonding started from the simply 
supported end) . 
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Figure 14. Instability domains glass-epoxy with free and 
simply supported ends (Dedonding started from the free end) . 
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Figure 15. Instability domains boron -aluminum with free and 
simply supported ends (Dedonding started from the free end) . 
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Figure 16. Instability domains graphite -epoxy with free and 
simply supported ends (Dedonding started from the simply 
supported end) . 
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Figure 17. Instability domains of simply supported graphite 
epoxy (Dedonding started from the center) . 




Figure 18. Instability domains of simply supported Kevlar 
epoxy (Debonding started from the center) . 
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Figure 19. Instability domains of simply supported glass-epoxy 
(Dedonding started from the center) . 




Figure 20. Instability domains of simply supported boron - 
aluminum (Dedonding started from the center) . 
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Figure 21. Instability domains of simply supported graphite 
epoxy (Dedonding started from the ends) . 




Figure 22. Instability domains of simply supported Kevlar 
epoxy (Dedonding started from the ends) . 
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Figure 23. Instability domains of simply supported glass-epoxy 
(Dedonding started from the ends) . 




Figure 24. Instability domains of simply supported boron - 
aluminum (Dedonding started from the ends) . 
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Figure 25. Instability domains graphite -epoxy with free and 
simply supported ends (Dedonding started from the free end) . 
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Figure 26. Instability domains Kevlar-epoxy with free and 
simply supported ends (Dedonding started from the free end) . 
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Figure 27. Instability domains boron -aluminum with free and 
simply supported ends (Dedonding started from the simply 
supported end) . 
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IV. CONCLUSIONS 



An interfacial debonding between a fiber and a supporting 
matrix made a significant effect on the critical static 
buckling load of the fiber. An initial debonding of size 10 
percent of the fiber length was most critical for a simply- 
supported fiber regardless of its location. Three different 
fiber materials embedded in the same epoxy matrix yielded 
almost the same percentage reduction of the critical buckling 
load even if their absolute magnitudes of buckling loads were 
different. For a free-simply supported fiber, a local 
debonding located at the free end was most critical for the 
buckling load. Any size of debonding located at any other 
position did not influence the critical buckling load. In 
addition, a fiber breakage around the center of fiber along 
with a partial interfacial debonding reduced the buckling 
loads more significantly. 

Interfacial debonding effects on dynamic instability of 
embedded fibers were generally more critical than those on 
static buckling loads. Debonding not only lowered instability 
domains but also widened the domains. Debonding effects were 
more significant for a free-simply supported fiber than for a 
fiber with both ends simply supported. A local debonding of 
size of 10 to 20 percent of the fiber length was most 
influential in affecting the domain of dynamic instability. 
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Debonding effects were not sensitive to changes in the 
location of debonding or boundary conditions of a fiber. 
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